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ABSTRACT 

Context. Nuclear star clusters (NSCs) are ubiquitous at the centers of galaxies. They show mixed stellar populations and the spectra of many NSCs 
indicate recent events of star formation. However, it is impossible to resolve external NSCs in order to examine the relevant processes. The Milky 
Way NSC, on the other hand, is close enough to be resolved into its individual stars and presents therefore a unique template for NSCs in general. 
Aims. Young, massive stars have been found by systematic spectroscopic studies at projected distances R < 0.5 pc from the supermassive black 
hole, Sagittarius A* (Sgr A*). In recent years, increasing evidence has been found for the presence of young, massive stars also atR> 0.5 pc. Our 
goal in this work is a systematic search for young, massive star candidates throughout the entire region within R ~ 2.5 pc of the black hole. 
Methods. The main criterion for the photometric identification of young, massive early-type stars is the lack of CO-absorption in the spectra. We 
used narrow-band imaging with the near-infrared camera ISAAC at the ESO VLT under excellent seeing conditions to search for young, massive 
stars within ~2.5 pc of Sgr A*. 

Results. We have found 63 early-type star candidates aiR < 2.5 pc, with an estimated erroneous identification rate of only about 20 %. Considering 
their A'-band magnitudes and interstellar extinction, they are candidates for Wolf-Rayet stars, supergiants, or early 0-type stars. Of these, 31 stars 
are so far unknown young, massive star candidates, all of which lie al R > 0.5 pc. The surface number density profile of the young, massive star 
candidates can be well fit by a single power-law (oc R^'^), with T = 1.6 + 0.17 at < 2.5 pc, which is significantly steeper than that of the late-type 
giants that make up the bulk of the observable stars in the NSC. Intriguingly, this power-law is consistent with the power-law that describes the 
surface density of young, massive stars in the same brightness range atR< 0.5 pc. 

Conclusions. The finding of a significant number of newly identified early-type star candidates at the Galactic center suggests that young, massive 
stars can be found throughout the entire cluster which may require us to modify existing theories for star formation at the Galactic center. Follow-up 
studies are needed to improve the existing data and lay the foundations for a unified theory of star formation in the Milky Way's NSC. 
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1. Introduction 

Nuclear star clusters (NSCs) are ubiquitous in galaxies and ap- 
pear as comp act clusters at the dynamical cen ters of th eir host 
galax ies (e.g.. iBoker eTa ni2002l: ICaroUo et al.l [l998: Co te et all 
l2006l) . They have luminosities in the range 10^-10*^ Lq, effective 
radii of a few pc, and masses of the order 10^ - 10*^ Mq. They 
are typically 1-2 orders of magnitude brighte r and more mas- 
sive than globular clusters (IWalcher et al. ' 2005), which places 
NS Cs among the most massive known clusters in the Universe 
(see lB6kerll2010l for a brief review on NSCs). Star formation 
in NSCs appears to be a (quasi-)continuous process. The major- 
ity of NSCs have mixed old and young stellar populations and 
show frequently signs o f star formation within the past lOOMyr 
(e.g.. IWalcher et akl 12006) . NSCs show complex morphologies 
and frequently coexist with supermassive black holes (SMBHs; 
ISethet al.ll2006ll2008l) . 

Sagittarius A* (Sgr A*), the SMBH at the center of the Milky 
Way (MW), with a mass of roughly 4 mill ion solar masses 
and located at a distance of about 8kpc (e.g.. lGhez et al.ll2008t 
[Oillessen et al. 2009), is surrounded by a dense and massive 
star cluster, the MW's NSC, which has an estimated mass of 
3 + 1 .5 X 10^ Mq ( "Launhardt et al.!l2002l) and a half light radius 
of 3-5 pc jGraham & Spider 2009; Sch6deil l20TTl see the latter 
paper for a brief review on the MW NSC). The mid-infrared im- 
ages from the Spitzer Space Telescope show how t he NSC stands 
out as a separate structure at the center of the MW dStolovv et al.l 
120061) . 



Star formation events occurred in the central parsec of the 
Galactic center (GC) about 10^ and a few times 10^ years 
ago (Krabbe et al. 1995). The massive, young stars formed in 
the most recent star formation event were found to be con- 
centrated within a projected radius of R - 0.5 pc around the 
central SMBH, with their density increasing toward Sgr A*. 
About half of them appear to be located within a disk-like struc- 
ture (Levin & Beloborodov 2003; Paumard et al. 2006; Lu et al] 

2009) . The currently favored hypothesis for their formation 
is that they formed about 6 Myr ago in situ in a dense gas 
disk around the SMBH (e.g.. lBonnell & Ricell2008l;lBartko et akl 

2010) . An alternative scenario for the presence of young, mas- 
sive stars in the central parsec is the formation of a massive 
cluster at more than several parsecs distance from the SMBH, 
followed by its infall toward the central parsec via dynami- 
cal friction and its subsequent dissolution (e.g., Gerhard 200ll 
Kim & Morris 2003). However, this hypothesis encounters dif- 
ficulties in explaining observations such as the steep increase 
of the surface density of the young stars toward Sgr A* and 
the relatively low total mass of the young, massive stars in 
the central 0.5 pc (see discussion in Bartko et al. 2010). Also, 
the apparent absence of X-ray radiation from young, low-mass 
stars in the central parsecs dis favors the cluster infall hypothesis 
(iNavakshin & Sunvaevll2005h . 

Extinction toward the central parsecs of the GC is extreme. 
Ay = 30 - 50 mag and variable on arcsecond scales (e.g., 
IScoville et al.ll2003l;[Sch6del et al.l20T^ . which means that stel- 
lar colors are dominated by extinction, and broad-band col- 
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ors can hardly be used to distinguish between late- and early- 
type stars (see, however, the carefully extinction coiTected color- 
magnitude diagram in Fig. 8 of Schodel et al. 2010). Alternative 
methods to identify young, massive stars are spectroscopy. X-ray 
plus infrared multi-wavelength observations, search for emis- 
sion lines, or near-infrared imaging with filters sensitive to the 
CO-absorption of late-type stars. Since spectroscopic observa- 
tions are very intensive in observing time, they were almost 
exclusively focused toward a region within R ~ 0.5 pc from 
Sgr A *, with only a few observed fields at 0.5 < R < 0.8 pc 
(e.g.. lDo etaf] 120091: iBartko et allbolOh . Despite of the small 
area examined, spectroscopic observations dominate our current 
understanding of star formation at the GC with the view that 
young, massive stars can almost exclusively been found within 
R ~ 0.5 pc of Sgr A*. 

However, evidence for the presence of young, massive stars 
in the NSC at distances R > .5 pc from Sgr A* h as been steadily 
increasing in the past years. iGeballe et alJ (12006) analyzed the 
spectrum of a single star, IRS 8, located about 1 pc north of 
Sgr A*, and identified it as a 05-06 giant or supergiant. This 
source was an obvious target for dedicated spectroscopic ob- 
servations because it is su rrounded by a bow-shoc k and thus 
stands out from the field. iMauerhan et alj (1201 Obi ) combined 
X-ray and near-infrared observations to identify a few dozen 
Wolf-Rayet stars and O giants or supergiants throughout the 
central region (at R < lOOpc) of the MW. These observations 
are probably particularly sensitive to massive bin aries, where 
the X-ray emission is produced in colliding winds. Wa ng et al] 
(|2010). Dongetal. (2011), and Mauerhan et al. (2010iJ) used 
Pa a narrow -band imaging with HST to find young, massive 
stars throughout within a few tens of pc of Sgr A*. This method 
picks out most likely massive stars with strong winds, but it runs 
into problems in the central parsecs near Sgr A* because of the 
bright emis sion from the interstellar medium in the so-called 
mini-spiral dEkers et al.lll983l:lLo & Clausseni ri983) and the rel- 
atively low angular resolution that they used (FWHM w 0.2"). 
An efficient way of distinguishing between young, massive stars 
and old, late-type stars in imaging data is near-infrared imaging 
that makes use of the CO-a bsorption at A ^ 2 . 3 um in the at- 
mosp heres of the cool stars dGenzel et alj|2003l: iBuchholz et alj 
l2009l) . Young, massive stars are identified in this method through 
the ab sence of CO -absorption combined with their high lumi- 
nosity. |Buc^olzetaL| (2009) used this method, combined with 
adaptive optics (AO) imaging and found evidence for a popu- 
lation of young, massive stars at R > 0.5 pc from the central 
SMBH. However, the use of AO with the coiTespondingly small 
instrumental field-of-view (FOV) restricted them to a relatively 
small region atR< 20"/0.8 pc from Sgr A*. 

Here, we present the results of a systematic search for young, 
massive stars within R » 2.5 pc of Sgr A* via imaging data from 
the VLT near-infrared camera ISAACQ. Our method is naiTow- 
band seeing-limited imaging under good seeing conditions, us- 
ing the CO-band absorption of late-type stars to distinguish be- 
tween young, massive, early-type stars and late-type giants in 
the central ~6 x 6 pc. Recurrence to seeing-limited observations 
allow us to probe a significantly larger FOV than what would 
be possible in comparable time with AO imaging. Our work 
has some overlap with, but is also complementary to the above- 
mentioned work. 



' Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile under pro- 
gramme 081. B-0099. 



2. Observations and Data Reduction 

Imaging observations were acquired through various narrow- 
band filters with ISAAC at the ESO VLT dMoorwood et al.l 
1998). The FOV is 2:5 x 2:5 with a pixel scale of a' 15. Details 
of the observations are listed in Table[T] The angular resolu- 
tion of observations with ISAAC are limited by seeing, which 
is of special importance here because the GC is an extremely 
crowded field. Note that the DIMM seeing measurements listed 
in Tab.[T]present a pessimistic over-estimation of the actual see- 
ing values. On the one hand, according to the Roddier formula, 
the seeing FWHMoc /l^*'^, hence visual seeing of 1.0" trans- 
lates to about 0.8" at 2.09 jum. On the other hand, the DIMM 
measurements on Paranal are influenced by a turbulent layer 
close to the ground that hardly affects the telescopes themselves 
dSarazin et alj2008l) . The average point-source FWHM in aU our 
images is therefore only about 0.5" - 0.6". 

Most of the data reduction process was standard, with sky 
subtraction, bad pixel correction, and flat-fielding. The sky ob- 
servations were acquired with identical settings as the target ob- 
servations, but on a dark cloud about 713" west and 400" north 
of Sgr A*. 

A non-standard part of the image reduction process was 
taking care of detector row bias variations as described in the 
Instrument Data Reduction Cookbook, issue 80.0, edited by ESO 
Paranal Science Operations and available on the ESO web site. 
Since different rows of the detector are read out at different times 
when applying double coiTelated read, different parts of the de- 
tector sample different parts of the non-linear integration ramp. 
This can lead to important biases for short integration times and 
bright sources. The bias variations depend on integration time 
and detector illumination, and are not constant in time. In par- 
ticular, the bias depends on the brightness and distribution of 
sources over the FOV and is therefore different for each point- 
ing. The bias variations are non-uniform across the detector and 
cannot be modeled by a smoothly varying function. Since the 
upper and lower halves of the detector are read out in parallel 
there is a visible jump of the background around detector row 
512 in individual exposures. 

However, the bias variations are largely restricted to the ver- 
tical (y-)direction on the detector and the bias is, to 0-th order, 
uniform along each row. Therefore, the bias can, in principle, be 
reliably determined by measuring the median along each row, 
after suitably clipping the highest and lowest pixel values. A dif- 
ficulty in this procedure is that the GC is extremely crowded and 
there are hardly any pixels on the detector that do not contain 
part of the PSF of a star. After some experimenting we chose to 
use the median of the 200 lowest pixels from each row to esti- 
mate the bias. The 200 lowest pixels usually fall on dark clouds 
in almost each row in each image, where the bias can then be de- 
termined with reasonable accuracy, but at the price of losing in- 
formation on the absolute background emission. After bias cor- 
rection the images were combined into a final mosaic. 

The 2.25 jum naiTow-band image is shown in Fig.[Tl along 
with zooms onto a detail in the final mosaics composed from the 
bias-corrected and from the uncorrected images. The procedure 
for bias correction worked well, but some systematic errors re- 
mained, particularly because the brightness and density of stars 
is highly non-uniform across the images, which leads to over- 
estimation of the bias in rows crossing the center of the cluster 
and under-estimation of the bias in rows that cross dark patches 
of strong extinction. Since the final mosaics were composed of 
dithered images, the systematic errors of the bias correction were 
somewhat smeared across the final mosaics. Since we used PSF 
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Table 1. Summary of observations. 



UTC Date 


Filter ID 




A/l* 


DIT" 


NDIT^' 


N" 


Seeing' 






yum 


yum 


s 






arcsec 


2008 May 19 


NB1.71 


1.71 


0.026 


5.0 


15 


5 


0.6- 1.1 


2008 Sep 17 


NB1.71 


1.71 


0.026 


5.0 


15 


5 


0.9 


2008 Jul 18 


NB2.09 


2.09 


0.02 


5.0 


14 


5 


0.6 - 0.8 


2008 May 25 


NB2.25 


2.25 


0.03 


5.0 


15 


5 


0.9- 1.1 


2008 Sep 17 


NB2.25 


2.25 


0.03 


5.0 


15 


5 


0.7 


2008 May 25 


NB2.34 


2.34 


0.03 


5.0 


15 


5 


0.8 


2008 Sep 17 


NB2.34 


2.34 


0.03 


5.0 


15 


5 


0.9 



Notes. Central filter wavelength. Filter width. Detector integration time. Number of DITs averaged per dither position by readout 
electronics. Number of dither positions. According to DIMM measurements. 



fitting for photometry, with a simultaneous estimation of the 
background, the impact of a variable background on photom- 
etry is somewhat mitigated. Nevertheless, as discussed below, 
we could observe color trends in our data, mainly along the y- 
direction, with a jump near the middle of the images. We believe 
that these color trends are related to the uncertainties of the bias 
correction. A near-infrared three-color (1.19//m, 1.71 //m, and 
2.25 jum) composite image of the central ~ 150" region is shown 
in Fig.m 



3. Data Analysis 

3.1. Photometry 

We performed source detection and point-spread function 
(PSF) fitting with the IRAF (Imagi ng Reduction and Analysis 
Facility fl and DAOPHOT packages ( IStetsonll987h . We used the 
daofind task to identify point sources, and the sources were then 
input for PSF fitting photometry to the allstar task. The PSF 
model was constructed using a "Pennyl" functiorQ- From 50 to 
100 sources were used to construct the variable PSF for each im- 
age. The PSF model was linearly variable over the image with 
terms proportional to pixel offsets in the x- and y-directions. 

The photometric calibration was carried out usi ng a photo- 
metric catalog of point sources (ISchodel et al.ll2010l) . In the cat- 
alog, H-, Ks-, and L'-band magnitudes of point sources in the 
central 40" x 40" region are listed. The PSF-fitting magnitudes 
for the NB1.71 filter were calibrated with the //-band magni- 
tudes, and those in the NB2.09, NB2.25, NB2.34 bands were 
calibrated with the Ks magnitudes. The small bias introduced by 
this procedure among the NB2.09, NB2.25, and NB2.34 bands 
could be neglected for our purpose because we are primarily 
interested in differential photometry, where overall offsets are 
of little importance. In the following, we search for early-type 
stars by using the magnitude difference between the NB2.25 and 
NB2.34 bands. 

An astrometric calibr ation was done us i ng the IRSF/SIRIUS 
point source catalog by iNishivama et al.l (l2006h . Positions of 
point sources in this catalog were calibrated with the 2MASS 
point source catalog (Skrutskie et al. 2006), and the positional 
accuracy of the IRSF/SIRIUS catalog was estimated to be ~ 
0.1". Due to image distortion, differences between stellar po- 
sitions in the IRSF/SIRIUS catalog and ISAAC images are as 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

^ A Gaussian core with Lorentzian wings. 



large as 0.3" at the corners of the images, although the differ- 
ences at the image center are within 0.1". This astrometric cali- 
bration was done in order for matching of our sources with other 
observations. The matches in stellar positions were calculated 
automatically at first, but a confirmation of the matching was 
done by eye, by plotting the position of the matched sources on 
the ISAAC images. 

To avoid spurious detections of point-like sources, the 
sources were used in the following analysis only if a source was 
detected in all of the five observations of the NB2.09, NB2.25, 
and NB2.34 bands (see Table[TJ. 

3.2. Color correction and source selection 

Although PSF fitting was done with optimized parameters and 
a variable PSF, the resulting color magnitude diagram still had 
a large scatter (see the left panel in Fig. [3]). The typical error in 
the [2.25] - [2.34] coloflis from 0.15 to 0.19 at [2.25] < 13.5, 
which is larger than the expected difference of ~ 0.1 mag 
between early-type d warfs and late-type giants at ~ 2.3/vm 
(iBuchhoIz et al.ll2009l) . To investigate the cause of the large scat- 
ter, we made [2.25] - [2.34] color distribution profiles as func- 
tions of the image x- and y- axes (upper panels in Fig. |4]i. We 
can see a clear color trend along the y-axis and a discontinuity 
at y 500 pixel, the probable cause of which is described in the 
last paragraph of ^ The mean [2.25] - [2.34] color at y ~ is 
0.1, but at y ~ 950 it is about H-0.4. This trend is thus identi- 
fied as the principle cause of the large scatter in the CMD. There 
is also a trend along the x-axis, but it is continuous and has a 
smaller amplitude. 

To improve the CMD, we carried out a color correction to 
reduce the influence of the color trend, and eliminated stars with 
large photometric uncertainties. For the color correction, we as- 
sumed that the average of the intrinsic stellar colors is the same 
throughout the observed field. This assumption is appropriate 
for our observed field due to similar intrinsic colors of almost 
all stellar types observable around the K band, the dominance of 
late-type giants in the GC region, and the restricted wavelength 
range of our observations. This assumption is also well justified 
in the case of GC, and ha s been used in previous work, e.g., 
Schodel et al.l (I20()7ll2010h . where this point is further discussed 
and illustrated. 

At first, to correct the color trend along the y-axis, the ob- 
served field was divided into sub-fields of 1024 (x) x 20 (y) pix- 
els. The [2.25] - [2.34] color histograms of stars in the sub-fields 
were constructed and fit with Gaussian functions. As described 



[A] denotes a magnitude in a narrow-band filter with the central 
wavelength of A. 
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above, the intrinsic colors of the histograms' peaks are expected 
to be almost the same, and the difference between the central 
wavelengths of the NB2.25 and NB2.34 filters is only 0.09 jum. 
Preliminary color-shifts for each sub-field were therefore calcu- 
lated by requiring that the peak of each Gaussian had a color of 
[2.25] - [2.34] = 0. After the color-shifts along the y-axis we 
carried out a similar procedure along the x-axis, by dividing the 
observed field into sub-fields of 20 (x) x 1024 (y) pixels. The 
[2.25] - [2.34] color distribution profiles after the color correc- 
tion are shown in Fig. |4l bottom panels. After correction there 
appears to be no significant color trend along either the x- or the 
y- axis. The CMD after the color correction is shown in the right 
panel in Fig. [3] 

Next, stars with large photometric uncertainties were elim- 
inated from our source list. The stellar density is extreme in 
the GC, which means that a large fraction of sources will in- 
deed be affected by systematic photometric errors due to source 
confusion. Comparing two images taken under different seeing 
conditions is an efficient way of identifying and excluding those 
sources from our sample. Since two independent sets of observa- 
tions existed for the NB2.25 and NB2.34 bands (see Table[T]), we 
could estimate the photometric uncertainties of the stars from the 
differences between these two measurements. We calculated the 
RMS of the photometric uncertainties in bins of one magnitude 
width, and eliminated stars with uncertainties larger than the 
RMS. Through this procedure, 3063 sources out of 8014 were 
excluded from our source list. The final CMD with 495 1 sources 
is shown in Fig. |5] The typical combined uncertainties of the 
[2.25] - [2.34] color (squared sum of the photometric uncertain- 
ties in both bands) now range from 0.06 to 0.08 at [2.25] < 13.5, 
more than a factor of two smaller than those without the color 
correction and the elimination of the sources with large photo- 
metric uncertainties. 

3.3. Source Classification 

The final CMD has sufficient photometric accuracy to differ- 
entiate early-type stars from late-type giants. In Fig. |5j we 
show the [2.25] - [2.34] CM D, in which sp ectroscopically con- 
firmed late-type giants (Ma ness et al.ll2007h and early-type stars 
dPaumard et alj 12006. Table B.l) are marked by green circles 
and red squares, respectively. The majority of sources in our fi- 
nal sample lie on a sequence from ([2.25] - [2.34], [2.25]) ~ 
(-0.2,9.5) to (0.1,15.5). This coincides with the location of 
most of the spectroscopically identified giants, indicating that 
the sequence is the red giant branch (RGB) on this CMD. On 
the other hand, at [2.25] ^ 12.5, most of the known early-type 
stars (green circles) are located to the right of the RGB, and only 
two of a total of 23 early-type stars lie on the RGB. Red color 
in [2.25] - [2.34] means no significant CO absorption at 2.34 
yum, thus suggesting an early spectral type. Hence Fig.|5]reveals 
that early-type stars can be identified reliably in our sample at 
[2.25] ^ 12.5. From the location of the spectroscopically iden- 
tified stars, we estimate a probability of false identifications on 
the order of 10 % in this region of the CMD. 

To automatically identify early-type stars and to determine 
the magnitude limit of our procedure, we proceeded as follows. 
We constructed [2.25] - [2.34] color histograms for stars within 
one magnitude wide bins in [2.25], and fit the histograms with 
Gaussian functions. The mean and standard deviation, cr, of the 
best fit for each bin are shown by pink circles and pink bars in 
Fig. |6] respectively. The mean and cr correspond to the mean 
color and width of the RGB in each magnitude bin. We defined 
stars more than 2 cr redder than (i.e., to the right of) the mean 



RGB color as early-type star candidates (light green crosses in 
Fig.©. 

The color histograms of the RGB are expected to follow a 
normal distribution. When we construct a histogram of a "rela- 
tive" [2.25] - [2.34] color from the mean of the RGB, the his- 
togram shows a normal distribution of the RGB plus an excess 
component of the early-type stars at the red side of the histogram 
(Fig-EJ- So the "blue" stars, which are the blue tail of the normal 
distribution of the RGB, can be used to estimate the contamina- 
tion of the early-type star candidates by RGB stars because of 
measurement errors. So we define stars more than 2 cr bluer than 
the RGB as "blue" stars (blue crosses in Fig.|6]l. 

Magnitude histograms of the early-type candidates and the 
"blue" stars are shown in Fig. [8] top panel. The number of 
blue stars is very small at [2.25] < 12.25, but it increases at 
[2.25] > 12.25. The bottom panel in Fig.[8]shows the ratio of the 
blue stars to the early-type candidates, i.e., the expected contam- 
ination of the early-type candidates by erroneous identification. 
The ratios are less than ~25% at [2.25] < 12.25, and more than 
60 % at [2.25] > 12.25. Therefore we can conclude that the 
magnitude limit above which we can differentiate the early-type 
star candidates from the red giants is [2.25] * 12.25. The stars 
with [2.25] < 9.75 are also excluded from the candidates due to 
unreliable photometry caused by detector saturation. 

In the magnitude range of 9.75 < [2.25] < 12.25, we 
have found 63 early-type star candidates. In the same magni- 
tude range, 20 spectroscopically confirmed early-type stars are 
included in the final CMD (green circles in Fig.|5]and|6l from 
Schodel et al.„2009„ Table B. 1). Of these, only one source is mis- 
identified as a red giant, and the others are recognized as early- 
type stars in our analysis. 

In the same magnitude range, we have found 12 "blue" stars 
located at the left side of the RGB. Assuming that the "blue" 
stars are red giants with a large photometric uncertainty, and the 
number of falsely identified early-type candidates is similar, the 
probability of contamination is estimated to be 12/63 = 0.19. 
However, as described below, 5 of the 63 candidates are not 
early-type stars, but intrinsically very red sources. By assuming 
that these two sources of uncertainty are independent, this in- 
creases the estimated contamination of early-type candidates to 
21 %. Note that this is a conservative estimate because the very 
red sources could be excluded by an additional color cut. This 
would, however, complicate our analysis unnecessarily. Also, 
some of these sources are probably early-type stars that are sur- 
rounded by bow-shocks (.Tanner et al...2005.) . 

3.4. [2.25] vs [1.71] - [2.25] CMD 

We constructed an extinction-corrected [2.25] vs [1.71] - [2.25] 
CMD (Fig. |9]l which is expected to be very similar to a broad 
band H vs H - K CMD. For the data sets of the NB1.71 and 
NB2.25 filter, we did the same procedure of the color correction 
and the elimination of the sources with large photometric un- 
certainties as described in ^3.21 Then, we carried out "relative" 
extinction correction following the procedure in ISchodel et al.l 
(2010); we calculated the mean of the [1.71] - [2.25] colors of 
the 10 nearest stars, ([1.71] - [2.25])mean, at each position, and 
the mean color excess -£'[i.7i]-[2.25] = ([1-71] - [2.25])niean - 2.1 
was used for the "relative" extinction correction of the sources 
in the NB1.71 and NB2.25 bands. Here we assume that most 
of the nearest stars are late-type red giants which have a mean 
color of [1.71] - [2.25] = 2.1 at the GC. More than 1 -magnitude 
bluer ([1.71] - [2.25] < 1.1) or redder ([1.71] - [2.25] > 3.1) 
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sources are excluded from the nearest star list because they are 
considered background or foreground stars. 

By comparing both the [2.25] vs [2.25] -[2.34] CMD (Fig.© 
and the [2.25] vs [1.71] - [2.25] CMD (Fig.©, a major benefit 
of the narrow-band observations with the NB2.25 and NB2.34 
filters is easily understandable. In the [2.25] vs [1.71] - [2.25] 
CMD, it is difficult to find a difference in the distributions be- 
tween t he spectroscopically confirmed ea rly-typ e stars (green 
circles; iPaumard et a J 20061: Schodel et al.l l2009l) and late-type 



giants (red triangles; iManess et al.ii2007i) even after our extinc- 
tion correction. This means that the accuracies of the photom- 
etry and extinction correction are not sufficient to distinguish 
the early-type stars from the giants. On the other hand, al- 
most all of the known early-type stars are identified as early- 
type star candidates in the [2.25] vs [2.25] - [2.34] CMD at 
9.75 < [2.25] < 12.25 (Fig.©. 

Although the [2.25] vs [1.71] - [2.25] CMD is not practical 
for the early-type star search, it is very useful to find foreground 
objects because the [1.71] - [2.25] color is sensitive to the in- 
terstellar extinction. Of the early-type star candidates, only one 
source has a very blue color of [1.71] - [2.25] « 0.7, and it 
is clearly a foreground object. However, the other sources have 
similarly red colors as the spectroscopically confirmed early- 
type stars and giants in the GC. They are thus likely to be lo- 
cated in the GC. To eliminate foreground contamination, stars 
with [1.71] - [2.25] < 1.5 are excluded from the following anal- 
ysis. 



4. Results and Discussion 

According to the definition of the early-type star candidates de- 
scribed above, we have found 63 candidates in the range of 
9.75 < [2.25] < 12.25 (Table Ell. Considering their magnitude 
and the interstellar extinction of Ak ~ 2 - 3 mag in this region, 
they are candidates for Wolf-Rayet stars, supergiants, or early 
O-type stars. 

An advantage of the early-type star search with the NB2.25 
and NB2.34 narrow-band filters is an almost negligible influence 
of interstellar extinction. The small difference (0.09jum) of the 
central wavelengths of the two filters introduces a very small 
color excess of £'[2.25]-[2.34] = 0.07 Aks, according to the inter- 
stellar extinction law toward the GC (A^ oc A" and a ~ -2.0, 
where A^ is the amount of extinction a t the wavelength A; 
iNishivama et alJl20 06': 'Schodel et ani2010h . In addition, the red- 
dening vector on the [2.25] vs [2.25] - [2.34] CMD is almost 
parallel to the RGB (Fig. ©. So the objects on the CMD show 
shifts parallel to the RGB due to the interstellar extinction, which 
has thus a negligible effect on the classification of the sources. 



4.1. Comparison with otiier observations 

Intrinsically very red sources which we here call dust embedded 
sources (DBS), exist i n the central parse c of our Galaxy (e.g . , 
iMoultakaetatlllooi iTanneret al.ll2005t iGebaUe et all l200a) . 
They are included in our list of the early-type star candidates 
because there is good evidence that at least some of them are 
young, massive stars ( Tanner et al. 2005). The position of the 
DESs on our CMD is shown in Fig.|6] With the exception of IRS 
2L, all of them lie far to the right in the CMD. Several known 
DESs, such as IRS IW, 9, 13, are excluded from our final source 
list because of insufficient photometric accuracy. 

In the region we observed there are 18 sources with sig- 
nificant Pa Of emission (the primary Pa a emitting sources in 



iDong et aDl201 ih . Since most of them are concentrated within 
0.5 pc of Sgr A*, where source confusion is severe in our im- 
ages, only seven sources are found to be matched wi th our final 
list (source ID 117, 121, 123, 124, 129, 131, 133 in lDong et all 
2011)0 The sources (122, 127, 128) are not identified as point 
sources due to source blending, one (128) is too faint for reliable 
detection, and the others (38, 118, 119, 120, 125, 126, 132) are 
excluded from our list due to large photometric uncertainties. As 
shown in Fig. |6] (purple squares), all of the seven sources in our 
list, one of which is as faint as [2.25] ~ 14, are distributed on 
the right side of the RGB. It follows that our method can reli- 
ably identify young, massive stars that were found via their Pa a 
emission. 

After excluding the DESs, we find that 20 of our early-type 
candidates lie within the fields observed with VLT/SINFONI 
(Bartko et al. 20100. Of those, 18 sources are confirmed as 
early-type stars with spectroscopic observations. This fact sup- 
ports our claim that the sources distributed on the right side of 
the RGB on the CMD are reliable candidates for early-type stars. 

Spectra of six of our early-type candidates were obtained 
by Blum et al. (2003), who identified four of them as hot/young 
stars (IRS 3, 6E, 21,15NE) and two of them as late-type giants 
(see Table|2]i. 

Twenty-eight of our early-type star candidates, after hav- 
ing excluded the DESs, have been previously observed with 
VLT/NACO and narrow-band filters (Buchholzet al.. .2009h . 
Although six of those sources were classified as late-type 
stars, 18 of them were listed as high quality early-type can- 
did ates. Two of th e remaining sources could not be classified 
by Buch holz etaP (2009) because of a noisy (IRS 16SW) or 
strongly reddened SED (ID 480), while a further two sources 
are not listed in their catalog. 

To summarize the comparison with past observations: we 
have found seven Pa a sources (Dong et al. 2011) in our source 
list, and all of them are identified as early-type star candidates; 
20 of our c andidates are included in the observed fields by 
iBartko et al ](I2010), and 18 o f them are spectros copically classi- 
fied as young, massive stars: Blum et al.| ( |2003l) observed spec- 
troscopically six of our candidates and report four of them as 
early-type stars; 18 among 28 matches between our candidates 
and st ars from the AO imaging observations bv iBuchholz et al.l 
(l2009l) are classified as high-quality early-type star candidates 
in their work, and only six of the matches are classified as late- 
type stars. The complete list of young, massive star candidates 
found in our observations, along with their measured colors, is 
contained in Table [J] 

4.2. Spatial distribution 

In Fig. [To] we show the NB2.25 filter image on which we have 
marked the early-type star candidates (light green circles) and 
the stars on the RGB (red dots) with 9.75 < [2.25] < 12.25. 
There is a significant concentration of the early-type candi- 
dates within a projected distance ^sgrA* of 0.5 pc from Sgr A*. 
These young, massive stars are well know n from previous imag- 
ing and spectroscopic observations (e.g.. iPaumard et all 120061; 
iDo et al.1 120091: iBartko et all l2010h . However, more flian half 
(35) of our candidates are found at i^sgrA* > 0.5 pc, in regions 



' # 121 is too bright, and # 131 and 133 are too faint for reliable 
identification of their spectral type, so they are not included in our list 
of early-type star candidates 

* Since their source list has not been published, we made a matching 
by eye with their Fig. 1. 
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that have so far hardly been examined with spectroscopic ob- 
servations. Although spectroscopic confirmation is necessary, 
the large number of the candidates, the small contamination 
probability of ~ 20 %, and the findings by other authors (e.g., 
iDong et al.ll20Tll) suggest the existence of a significant number 
of young, massive stars at distances R > 0.5 pc from Sgr A*. 

The intrinsic distribution of young, massive stars will of 
course differ from the observed one because the probability of 
identifying a candidate at a given location depends on factors 
such as extinction, crowding and local PSF. The apparent lack 
of sources in the south eastern and north western parts of the 
FOV is due to ext remely high extinction in these regions (see 
ISchodel et al.ll2007i) . Due to strong PSF distortion near the cor- 
ners and edges of the FOV of our observations, the density of 
sources in our sample (both the red giants and early-type stars) 
decreases significantly toward the south-eastern, north-western, 
and north-eastern corners of the FOV. However, none of these 
difficulties affects our main finding of a significant number of 
so far unknown candidates for young, massive stars at projected 
distances between 0.5 and 3.2 pc from Sgr A*. 

Azimuthally averaged, projected stellar surface density plots 
as a function of the distance to Sgr A* for all stars, red giants, 
and early-type candidates with 9.75 < [2.25] < 12.25 are shown 
in Fig.[TT] This analysis is limited to the region within 60" from 
Sgr A*, to extract stellar surface densities in full rings, and to 
exclude the regions with low completeness near the corners of 
our images. The bin widths were chosen such that each bin con- 
tains a constant number of stars (five stars for the analysis of 
early-type star candidates only, and 15 stars otherwise) in order 
to give each data point roughly the same statistical weight. 

The surface density profiles show clearly that the spatial dis- 
tribution of young, early-type stars is different from that of old, 
red giants outside the central 0.5 pc. We find an apparently con- 
tinuous surface density profile for the early-type star candidates 
in the range r.'5 < ^sgiA* ^ 60". The profile was fit with a power 
law (oc R^^) which gives a spectral index of FEariy = 1 -60 + 0. 17. 
On the other hand, the profiles of the entire sample and for the 
red giants show a break around 10" (as known from previous 
work, e.g., lEckart et aUfToollSchodel et alJl2007h . so the fitting 
was caiTied out only in the range of 10" < Z^sgrA* ^ 60", provid- 
ing power-law indices of Fan = 0.92±0.11 andFRG - 0.76+0.10 
for all stars and red giants, respectively. For comparison, we 
also fit the profile for the early-type star candidates only for 
RsgvA* - 10" to 60", and obtain r^,^^.^^ ^jq// = 1.75 + 0.49. Note 
that we have not made any extinction and completeness correc- 
tions for these plots, and thus the absolute values of F may have 
additional systematic errors; however, it is possible to compare 
the three indices, T^n, Frg, and FEariy (or Fg^j.jy ^j^//), because 
they suffer almost the same interstellar extinction and have the 
same photometric completeness limits. Also, the power-law in- 
dices inferred here agree well with previous work (see below). 
The early-type star candidates show a steep, apparently continu- 
ous radial profile, while the red giants show a shallower profile 
with a break at ~ 10". The significant difference of the profiles 
also suggests that not every early-type star found in our analysis 
is contamination by erroneous identification of a late-type giant. 

While the power-law indices derived here may suffer from 
some additional systematic uncertainties as discussed above, 
those for the entire sample and the red giants, Fau - 0.92 and 
Frg - 0.76 a r e in g ood agreement with previous observations. 
ISchodel et al.l (l2007h produced a fit of a broken power law to 
the surface number density profile of the stars in the NSC, re- 
gardless of type, and found F = 0.75 ±0.10 outside a break ra- 



dius at 6"0 + 1"0. Th rough narrow-band imaging observations, 
iBuchholz et aP (l2009h obtained F = 0.86 ± 0.06 and 0.70 + 0.09 
for all sources and late-type stars, respectively, at 6" < /?sgiA* ^ 
20". Here our limiting ma gnitude, [2.251 - 12.25, is much 
shallower than K ^ 1 7.75 iSchodelet al.ll2007h and K = 15.5 
(iBuchholz et al.l l2009t) . So the agreement implies that incom- 
pleteness and variable extinction do not affect our result signifi- 
cantly within the photometric limit of this work. 

Interestingly, the power-law indices for the early-type star 
candidates are also similar to previous results at /?sgrA* ~ 10". 
We construct a plane-of-the-sky projected surface density plot 
for early-type star s with 9.75 < K < 12.25 catalogued in 
iBartko et al.l (l2009h . and find Feaiiy = 1-47 + 0.42. This is con- 
sistent, within the uncertainties, with our values of ^^.^^ly^iQ" = 
1 .75 + 0.49, and also Fnai-iy = 1 .60 + 0.17. ^' 

4.3. Implications for Star Formation at the Galactic Center 

In the central parsec, where the potential is dominated by the 4 x 
IQ^ Mq SMBH, the densities for a molecular cloud to be stable 
against tidal disruption are on the order of a few 10'*^ (10**, 10^) 
nmcm'^ at r - 0.1 (0.5, 1.0)pc distance from SgrA*. These 
extreme densities are several orders of magnitude larger than 
what is fo und in m olecular clouds near the GC (see, e.g., dis- 
cussion in lGenzefe t al. 2003). Standard in situ star formation is 
therefore considered to be impossible in the central parsec of the 
MW. Two alternative scenarios were therefore suggested to ex- 
plain the presence of young, massive stars in the central parsec: 
(1) formation of the stars in a fra gmenting gas disc or streamer 
around the SMBH (e.g., iNavaks hin et al. 2007); (2) formation 
of a massive, dense stellar cluster at more than several parsecs 
distance and its subsequent infall and dissolution (e.g.. iGerhardl 
i2001h . 

The young, massive stars at the GC were found to be concen- 
trated mainly within R ~ 0.5 pc, with the highest concentration 
near Sgr A* and a steep fall-off of their surface density toward 
larger distances. Along with some e vidence for a very top-heavy 
initial mass function of these stars (iBartko et al . 2010, see how- 
ever Do et al., in preparation) and with the absence of X-ray 
emission from low-mass pre mai n sequence stars from a poten - 
tial dissolved progenitor cluster dNavakshin & Sunvaevll2005h . 
this appears to favour the scenario of star formation in a gaseous 
disk around SgrA*. 

It has to be kept in mind, however, that our current ideas 
on star formation at the GC are influenced by the limitations of 
observational studies. Those studies have been focused largely 
on the region within a projected radius of R < 1 pc from SgrA*. 
The main reason for this limitation is that spectroscopic searches 
are very costly and are therefore limited to small numbers of 
stars and/or to small areas. Therefore it is still unclear whether a 
clear outer edge of the distribution of the early-type stars exists 
or not. Recently, further spectroscopic observations foun d about 
10 ne w early-type stars outside the central 0.5 pc (Bartk o et al.l 
i20T0l) . and as described in the introduction, ever more observa- 
tional results support the existence of young, massive stars at 
larger distances from the SMBH. 

Here we present the results of the first systematic search, 
via narrow-band imaging, for young, massive stars within a pro- 
jected distance of R ~ 2.5 pc of SgrA*. Although our survey is 
limited to the brightest {K < 12.25) stars in the field because 
of the limited photometric accuracy imposed by the seeing- 
limited observations in this heavily crowded field, we detect a 
significant number of candidates for young, massive stars dis- 
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tributed throughout the studied region, with half of them located 
at R > 0.5 pc. The brightness range of these candidates sug- 
gests that they are of similar type - and therefore age - than 
the young, massive stars in the central 0.5 pc. As Fig. 10 shows, 
these stars appear to be randomly throughout the NSC (small 
number statistics keep us from studying this aspect more quan- 
titatively). Intriguingly, a single power-law, with an exponent of 
r - 1.60 + 0.17, can provide a good fit to the surface density of 
the young, massive stars as a function of distance from Sgr A* 
throughout the studied region. This power-law agrees within the 
uncertainties with the power-law that describes the surface den- 
sity of correspondingly bright young stars within R = 0.5 pc 
of Sgr A*, that w ere identified in spectroscopic studies (e.g., 
iBartko et al.ll2009l) . 

Our finding that a single power-law can describe the sur- 
face density of early type candidates throughout the examined 
region may indicate a physical connection between the young, 
massive stars inside and outside of the central 0.5 pc. Can these 
stars have formed in the outer parts of the hypothetical gas disk? 
It seems to be difficult to form massive stars via this mechanism 
at r > 0.5 pc. Hobbs & Navakshin (2009) suggest that the young 
stars at the GC could have been formed in a high-inclination 
collision of two massive gas clouds. In their simulation, gaseous 
streams or filaments, that extend as far out as 1 pc from Sgr A*, 
become self-gravitating and form stars. However, most of the 
stars formed in the streamers are low mass stars (0.1-1 Mq), 
while the most massive stars are concentrated towar d the central 
0.5 pc region, in contrast with our observations. Mapelli et al.l 
(120 12h also show that in the outer part (r > 1 pc) of the hypo- 
thetical disk, the gas density is lower and thus the Jeans mass be- 
comes about an order of magnitude larger than in the central few 
tenths of a parsec. Hence, the resultant star formation efficiency 
is very low outside a few tenths of a parsec. Since most of the ex- 
isting theoretical work has focused on explaining the presence of 
young, massive stars at r < 0.5 pc, without having observational 
evidence for a significant population of young, massive stars at 
larger distances, they were not optimized to explain star forma- 
tion at greater distances. Therefore, further studies are necessary 
to investigate whether star formation in infalling clouds/gaseous 
disks can explain the presence of young, massive stars beyond 
r ~ 1 pc as well. Also, deeper photometric studies, similar to the 
one described here, are needed to set tighter constraints on the 
statistics of young, massive stars throughout the region studied 
here and to confirm or reject the hypothesis that a single power- 
law can be applied to their surface density from to 2.5 pc. 

As concerns the cluster-infall scenario, our results do not sig- 
nificantly alter the currently accepted picture because the fun- 
damental observation stays valid that the young, massive stars 
are strongly concentrated toward Sgr A*. Nevertheless, we can- 
not yet reliably exclude the possibility that the candidates for 
young, massive stars at R > 1 pc may have been deposited at 
their current location by the dissolution of a stellar cluster at 
larger distances, in an ev ent that was unrel ated to the stars within 
R ~ 0.5 pc of Sgr A*. iGiirkan & Rasiol (12005), for example, 
showed that a ~ 5 x 10^ M0 stellar cluster with an initial position 
of lOpc from the SMBH can sink into the central a fewpc within 
a lifetime of young, massive stars (a few Myr). Since the surface 
density profiles predicted from the inspiralling cluster scenario 
strongly depen d on models employed in simulations (see, e.g. , 
simulation s by iKim & MorrisI 120031; iBerukoff & Hansenll2006t 
iFuiii etal.l l26l0r it is difficult to exclude the cluster scenario 
by comparing the surface density profiles. Kinematic informa- 
tion of the candidates is presumably crucial to discriminate the 
above scenarios. 



We can also consider the possibility that young, massive 
stars at R > 1 pc were formed in a different event (or events) 
than the ones at R < 0.5 pc. The required density of molec- 
ular clouds against tidal destruction by the tidal shear from 
Sgr A* is reduced to a few 10^nH2cm"^ at r = 2pc. Some 
clumps of molecular gas in the circumnuc lear disk (CND) 
may in fact reach the necessary densities (Christopher et alj 
2005; Montero-Ca stano et al. 2009) and some observations may 
indicate the onset of st ar formation in parts of the CND 
(lYusef-Zadeh et al.l l2008l) . It may thus be possible that the 
young, massive candidates formed at a distances from Sgr A* 
that were close to the currently observed ones. Recent surveys 
find in fact young, massive stars throughout the central tens of 
parsecs of the MW, indicating that massi ve stars may continually 
form throughout the central region (e.g.. lMauerhan et"al]|2010al: 
Wang et al. 2010). 

Currently, we do not have enough observational data to dis- 
criminate between the above scenarios. The direction of future 
work is clear. The immediate next step is to confirm (or dis- 
card) our candidates as early-type stars and determine their exact 
properties through spectroscopic observations. This will provide 
not only firm evidence for the existence of the young, massive 
stars outside the central 0.5 pc, but also mass and age informa- 
tion to understand their formation mechanism and their possi- 
ble relation with the young stars at R < 0.5 pc. Deeper imaging 
and spectroscopic observations will then be needed to search for 
fainter young, massive stars and establish their surface density as 
a function of location. Finally, proper motion measurements can 
provide further clues as to the nature and origin of the young, 
massive stars in the NSC. 



5. Conclusions 

We have carried out a systematic search for young, massive stars 
within a projected distance of R - 2.5 pc from Sgr A*, by means 
of near-infrared narrow-band imaging with VLT/ISAAC. Our 
method makes use of the absence of the CO-bandhead absorp- 
tion feature, which is ubiquitous in late-type stars. It is partic- 
ularly prominent in the red giants that almost completely dom- 
inate the brightness range of our stellar sample, which is lim- 
ited to the magnitude range of 9.75 < [2.25] < 12.25. We have 
found 63 early-type candidates at projected distances R < 2.5 pc 
from the Milky Way's central black hole. Through comparison 
with previous spectroscopic studies, which were mainly limited 
to R < 0.5 pc and via statistical considerations, based on the 
uncertainties of the measured stellar colors, we estimate the per- 
centage of erroneously identified stars as low as ~20%. About 
half of the candidates found in this study have not been previ- 
ously identified. All of those new candidates for young, massive 
stars are located atR > 0.5 pc. Our study thus shows that young, 
massive stars can be found throughout the nuclear star cluster. 
Follow-up studies, both spectroscopically and via AO assisted 
narrow-band imaging, will be needed to confirm in detail and 
deepen the findings of this work. So far, we lack a consistent 
model of star formation at the Galactic Center that can explain 
the presence of young, massive stars throughout the nuclear star 
cluster. 
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Table 2. Early-type star candidates. 
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Notes. Typical positional uncertainty is C.'l - Qf.'S (see 93.11 . Distance from Sgr A*. '''' Magnitude difference between [1.71] and [2.25], 
after the color correction ( 93. 2t . Magnitude difference b etween [2.25] and [2.34], after the color correction ( 93. 2t and the relative extinction 
correction ( 93.4t . '''' Name(s) of counterparts in lPaumard et al. (2006) and (Lu et al. 2009). The source # 53 is probably a foreground star, because 
of its very blue color in [1.71] - [2.25]. Spectral type from Paumard et al. (2006), while the superscri pt '(B a)' indicates the spectral type (early 
or no n-early) from Bartko et al. (201^, '(Ta) ' from Tanner et al. (2005), and '(Bl)' from Blum et al] i2003l) . Counterparts from Blum et aj] 
([2003, superscript 'Bl'),|Buchholz et al. (2009, superscript 'Bu') and Dons et al. i2011l . superscript 'Do'), and their source ID. El, E3, LI, L2, N, 
and R denote quality 1 early-type candidates, quality 3 early-typ e candidates, q u ality 1 late-type candidates, quality 3 late-type candidates, sources 
with noisy SED, and strongly reddened sources, respectively, in lBuchholz et al] ( l2009h . 
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m m m 

RA ['*] 



Fig. 1. Image of the GC through the 2.25 yum narrow-band filter of ISAAC (left). The inserts on the right show a detail after (top) 
and before (bottom) correction of the row bias, which removes most, but not all, of the variable detector background. 
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Fig. 3. [2.25] vs [2.25] - [2.34] color magnitude diagrams before (left) and after (right) the color correction. See g3.2l in detail. 
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Fig. 4. Color dependence of point sources along the x- and y-axes on the infrared array. We can see a clear color trend before the 
color correction (upper panels). However, after the correction (bottom panels), no clear dependence can be seen. See ^3.2l in detail. 
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Ear y type Stars 



Late type giants 
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Fig. 5. [2.25] vs [2.25] - [2.34] color magnitude diagram after the color correction and the elimination of sources with large 
photometric uncertainties. Known early-type stars and late-type giants are marked by green circles and red triangles, respectively. 
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Fig. 6. [2.25] vs [2.25] - [2.34] color magnitude diagram. Red color (positive value) in [2.25] - [2.34] means no significant CO 
absorption at 2.34 yum which, together with the brightness, is an indicator for massive, early-type stars. Typical uncertainties are 
represented by grey crosses at the left side. Stars with photometri c uncertainties large r than these typical ones are not included in 
this diagram. The reddening vector with a slope of A,\ oc A^^-^ (.Nishivama et alJl2006l) is shown at the upper right in the diagram. 
Dark green circles represent spectroscopically identified early-type stars. The means and sigmas of the RGB in the [2.25] - [2.34] 
color within bins of one magnitude width are shown by pink circles and bars, respectively. We define stars more than 2-cr redder 
than the red giants as early-type star candidates (light green crosses). Blue crosses, which are stars more than 2-cr bluer than the 
RGB, are used to estimate the false detection likelihood of the early-type star candidates. 
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Fig. 7. Histogram of stellar colors relative to the RGB mean (pink circles in Fig.|6l). Only stars with [2.25] < 12.25 mag are used 
for this histogram. A Gaussian fit is shown by a red line. An excess is clearly seen on the right (redder) side of the RGB, indicating 
the existence of early-type stars in this region. 
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Fig. 8. Top: Magnitude histograms for early-type star candidates (green) and stars more than 2-cr bluer than the RGB (see Fig.|6]l. 
This shows that the number of stars at the redder side (earlier spectral type) of the RGB is clearly larger than that at the bluer side 
(later spectral type). Bottom: Ratio of the bluer stars to the early-type star candidates. If we assume that the bluer stars are red giants 
with a large photometric uncertainty, a similar number of stars may contaminate our sample of the early-type star candidates. So 
this diagram shows the expected ratio of such contamination for the candidates. The uncertainties in the top panel are given by the 
square root of the counts in each bin. For stars brighter than [2.25] = 12.25 mag, the contamination is low (< 25%). 
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Fig. 9. [2.25] vs [1.71] - [2.25] color magnitude diagram. The early-type star candidates are overplotted by light green crosses. Dark 
green circles and red triangles represent spectroscopically identified early-type stars and late-type giants, respectively. Most of the 
candidates have a [ 1 .7 1 ] - [2.25] color (almost equivalent to aH -K color) of more than 1 .5, suggesting that they suffer from strong 
{Ak > 2) extinction, and thus they are not foreground sources. 
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Right ascension 

Fig. 10. Spatial distribution of the early-type star candidates with 9.75 < [2.25] < 12.25 (green circles) within the FOV of our 
observations. The candidates which have been unknown so far are indicated by blue arrows. Red giants identified by our analysis 
are marked by red circles. The large blue circle delimits a region within 0.5 pc (12'.'9) in projection from Sgr A*. There is a strong 
concentration of the early-type star candidates in the central 0.5 pc, as already known from previous studies. In addition, we have 
identified about 30 new candidates outside the 0.5 pc region. Note that the density of classified sources decreases strongly toward 
the north-eastern, north-western , and south-eastern corners due to unreliable photometry in these regions, caused by strong PSF 
variability. 
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Fig. 11. Plots of the azimuthally averaged, stellar surface number densities as a function of projected distance from Sgr A* for all 
stars in our sample (black), red giants (red), and early-type star candidates (green), with 9.75 < [2.25] < 12.25. The densities were 
determined in annuli of variable width, chosen such that each one contained a fixed number of stars (5 stars for the early-type star 
candidates, and 15 stars for the entire sample and the red giants). Data points at 10" < /?sgrA» ^ 60" were used for power-law fits 
of the radial profiles for our entire sample (black line), red giants (red line), and early-type star candidates (green line). The green 
dashed line represents a power-law fit to the densities of the early type candidates, including those at ^sgiA* ^ 10". 
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